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Abstract 

It is postulated that secreted soluble factors are important contributors of bystander effect and adaptive responses 
observed in low dose ionizing radiation. Using multidimensional liquid chromatography-mass spectrometry based 
proteomics, we quantified the changes of skin tissue secretome - the proteins secreted from a full thickness, reconstituted 
3-dimensional skin tissue model 48 hr after exposure to 3, 10 and 200 cGy of X-rays. Overall, 135 proteins showed statistical 
significant difference between the sham (0 cGy) and any of the irradiated groups (3, 10 or 200 cGy) on the basis of Dunnett 
adjusted t-test; among these, 97 proteins showed a trend of downregulation and 9 proteins showed a trend of upregulation 
with increasing radiation dose. In addition, there were 21 and 8 proteins observed to have irregular trends with the 10 cGy 
irradiated group either having the highest or the lowest level among all three radiated doses. Moreover, two proteins, 
carboxypeptidase E and ubiquitin carboxyl-terminal hydrolase isozyme LI were sensitive to ionizing radiation, but relatively 
independent of radiation dose. Conversely, proteasome activator complex subunit 2 protein appeared to be sensitive to the 
dose of radiation, as rapid upregulation of this protein was observed when radiation doses were increased from 3, to 10 or 
200 cGy. These results suggest that different mechanisms of action exist at the secretome level for low and high doses of 
ionizing radiation. 
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Introduction 

Both natural and man-made sources of ionizing radiation 
contribute to human exposure and consequently pose a possible 
risk to human health. While natural background radiation is 
unavoidable, the increased medical use of radiation, e.g. X-ray 
and CT scans for diagnosis, inevitably increases the public's health 
risk concerns over exposure of low doses of ionizing radiation. [1,2] 
Traditionally, ionizing radiation-induced cancer risks at low doses 
have been extrapolated from high dose studies, for which the 
effects on human health and biological systems are well 
documented, based on a linear no threshold dose-response model 
between radiation dose and effects; [3-5] however, recent literature 
suggests that the mechanisms of action for many biological 
Midpoints may be different at low doses from those observed at 
high doses. [6-1 1] 

Secreted proteins play key roles in cell-cell signalling, commu- 
nication, migration and growth, as they reflect the various stages of 
biological and pathological conditions. [12] Secreted soluble factors 
have long been postulated as the contributor of bystander effects 
observed from both in vitro and in vivo studies following ionizing 
radiation, where detrimental health effects of radiation exposure 



occur on neighbouring cells of the irradiated cells. [13] Many past 
studies on secreted factors were focused on a limited number of 
known proteins using targeted biochemical assays of Western-Blot, 
ELISA or protein microarrays;[14,15] despite decades of study, 
the identities of these secreted factors remain elusive. 

Using advanced LC-MS based proteomics technologies, we now 
present a comprehensive profiling of proteins secreted into 
culturing medium by a model 3D human skin tissue upon various 
doses of X-ray irradiation. Proteomic coverage of skin secretome 
was broadened by removal of medium proteins and extensive 
fractionation at peptide level before LC-MS/MS analysis. 
Quantitative comparison between low (3, 10 cGy) and high 
(200 cGy) doses of ionizing radiation indicates that different action 
mechanisms exist between low and high doses. 

Materials and Methods 

Tissue Culture and Irradiation 

The human skin tissues, EpiDerm Full Thickness 400, were 
purchased from the MatTek Corporation (Ashland, MA). Upon 
arrival, the tissues were randomized into 6-well plates, and 
cultured in EFT-400-MM maintenance medium (Dulbecco's 
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modified Eagle's medium supplemented with epidermal growth 
factor, insulin, hydrocortisone and other proprietary stimulators of 
epidermal differentiation) that also obtained from the MatTek 
Corporation. The tissues were grown in an atmosphere of 5% 
CO2 at 37°C and given 5 mL fresh medium every 2 days. 
Radiation experiments were performed after the tissues were 
acclimated for 5-7 days in our laboratory. 

The tissues were randomized into four radiation dose groups (0, 
3, 10 and 200 cGy) with twelve tissues (n= 12) per group. The 
irradiation of tissues with Low-LET (2 keV/ urn) X-ray was carried 
out using a Pantak XRAD 320 irradiator (Precision X-Ray, North 
Branford, CT) operated at 300 kV. The dose rates were as follows: 
3 and 10 cGy at 6 cGy/min, 200 cGy at 60 cGy/min. Following 
irradiation, tissues were returned to the same atmosphere as before 
collection at 37°C. At 48 hr post-irradiation, the media were 
collected and snap-frozen in liquid nitrogen and stored at — 80°C 
until further processing. 

Proteomic Sample Preparation 

Aliquots of each medium sample were pooled according to the 
radiation dose group, then 6 mL of the pooled samples and 
1.5 mL of each individual sample were subjected to immunoaffi- 
nity depletion to remove albumin and other abundant proteins 
using an IgY-14 LC5 column (Sigma- Aldrich) coupled with an 
Agilent 1100 series HPLC, under the conditions suggested by the 
manufacture for mobile phases and flow rate. The flow-through 
fractions (low abundance proteins) were collected and concentrat- 
ed in Amicon Ultra- 15 concentrators (Millipore) with MWCO of 
3 kDa, followed by a buffer exchange to 50 mil NH4HCO3 in the 
same unit according to the manufacturer's instructions. Sample 
proteins were next sequentially denatured with 8 M urea, reduced 
with dithiothreitol, alkylated with iodoacetamide, and digested 
with trypsin (Promega) at a trypsin/protein ratio of 1:50; the 
peptide mixtures were then cleaned with C 18 SPE cartridges 
(Sigma-Aldrich) and dried in vacuo before reconstituted with 0. 1 % 
formic acid for LC-MS analysis. 

Instrument Analysis 

For the pooled samples, protein digests were on-line fractioned 
and further separated using a custom-built 2D-LC system. The 
2D-LC system was composed of two Agilent 1200 nanoflow 
pumps and one 1200 capillary pump (Agilent Technologies), 
various Valco valves (Valco Instruments Co., Houston, TX) and a 
PAL autosampler (Leap Technologies, Carrboro, NC). Full 
automation was made possible by custom software that allows 
for parallel event coordination and therefore near 100% MS duty 
cycle through the use of one strong cation exchange (SCX) column 
(15 cmxl50|lm, 5|X PolySulfoethyl A particles, Poly LC Inc.) in 
the 1 st dimension, and two trapping columns (4 cmx lOOum, 3.6(1 
Aries C 18 particles, Phenomenex) and two analytical columns 
(35 cmx75u,m, 3n Jupiter C 18 particles, Phenomenex) in the 2 nd 
dimension. All columns were manufactured in-house by slurry 
packing of particles into fused silica (Polymicro Technologies Inc., 
Phoenix, AZ) using a 1 cm sol-gel frit for particle retention. Mobil 
phases used in the 1 st dimension were 0.05% acetonitrile in water 
and 500 mM ammonium formate to fractionate the samples 
(injection amount 9 u.g) into 15 fractions; and 0.1% formic acid in 
water and 0.1% formic acid in acetonitrile were used as mobile 
phases in the 2 nd dimension for reversed phase separation of the 15 
peptide fractions with a gradient length of 100 min. The eluents 
from the 2 nd dimension separation were electrospray ionized and 
analyzed by an Orbitrap Velos Pro hybrid mass spectrometer 
(ThermoScientific), which was operated in data-dependent MS/ 



MS mode, with one high resolution MS scan followed by 10 low 
resolution CID MS/MS scan events. 

For each individual sample, protein digests were separated on a 
4-column custom-built capillary LC system similar as reported 
previously. [16] Reversed-phase separations were carried out on 
35 cm x 75 um i.d. fused silica columns packed in house using 3|i 
Jupiter C 18 particles (Phenomenex). Mobile phases consisted of 
0.1% formic acid in water (A) and 0.1% formic acid in 100% 
acetonitrile (B) with a 100 min gradient. MS analyses were 
performed using an Exactive Orbitrap mass spectrometer 
(ThermoScientific) outfitted with a custom electrospray ionization 
(ESI) interface. MS spectra (AGC 3 x 10 6 ) were collected in 
positive ionization mode from 400 to 2000 ml z at a resolution of 
100 k. Triplicate analyses were performed on each individual 
sample. 

Data Analysis 

Identification and quantification of peptides was performed 
using the accurate mass and time (AMT) tag approach. [17] The 
datasets acquired in the 2D LC-MS/MS analyses were searched 
using the SEQJJEST algorithm to match the MS/MS fragment 
spectra with sequences from the UniProt human protein database 
(2010-05-05 version with 20,276 protein entries) with static 
carbamidomethylation of cysteine used as the modification. 
DeconMSN was used to accurately determine the parent ion 
monoisotopic mass. [18] Database-matched results were filtered 
using a very rigorous MSGF cut-off score of 1 E-12,[19] which 
resulted in an estimated false discovery rate of 0.3% at the unique 
peptide level in a decoy database search. The list of confidendy 
identified peptides was populated into an AMT tag database. The 
peptide elution times from each LC-MS/MS analysis were 
normalized to a range of 0 to 1 using a predictive normalized 
elution time (NET) model. [20] Both calculated monoisotopic 
masses and observed NETs of identified peptides were included in 
the AMT tag database. For the high resolution LC-MS datasets 
obtained from individual study samples, Decon2LS was used for 
peak finding and deisotoping at the individual MS spectrum 
level. [21] The output of Decon2LS that contains features of 
monoisotopic mass, intensity of most abundant isotope and elution 
time were further processed with VIPER to cluster mass and 
chromatographic elution features across datasets, [22] and to 
match the clustered features with the AMT tag database to 
identify the peptides and generate abundance profiles. Uniqueness 
probability (UP) score of >0.5 and an FDR of <10% were used 
for each dataset, and 1 ppm mass measurement error in at least 1 
dataset were used to remove the ambiguities in database matching 
and to improve the subsequent identification and quantification of 
clustered features. 

Crosstab file that contains peptide ID and peptide abundance in 
each sample was exported to MatLab® for statistical analysis. This 
includes quality control, normalization, protein roll-up, and 
peptide and protein level comparative statistical analyses. Peptide 
abundances were transformed to the logi 0 scale. Missing data 
values were not imputed. Quality control processing was 
performed to identify and remove contaminant proteins, redun- 
dant peptides, peptides with an insufficient amount of data across 
the set of samples, [23] and LC-MS runs that showed significant 
deviation from the standard behavior of all LC-MS analyses. [24] 
Peptides were normalized using a statistical procedure for the 
analysis of proteomic normalization strategies (SPANS) that 
identifies the peptide selection method and data scaling factor 
which introduces the least amount of bias into the dataset. [25] The 
peptide abundance values were normalized across the technical 
replicates with a peptide subset based on the L-order statistic using 
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a median centering of the data. Normalized logio abundance 
values of peptides were averaged across the technical replicates 
within each biological sample, and evaluated with a Dunnett 
adjusted 2-sided t-test to identify quantitative significance patterns 
in the peptide data. For protein quantification, peptide level 
significance patterns based on the comparisons of the 3, 10 and 
200 cGy groups to the control group (0 cGy) were used to select 
appropriate peptides for protein roll-up, which was based on the 
down-selection of peptides that have a common statistical trend 
followed by a standard average to a reference peptide, where the 
reference peptide was defined as the peptide with the least amount 
of missing data. Similarly, a Dunnett adjusted 2-sided t-test was 
used to assess differences in protein average abundance between 
the control group (0 cGy) and the 3, 10 and 200 cGy irradiated 
groups. Dunnett test p-values less than 0.05 were deemed 
statistically significant. The Uniprot data base (www.uniprot.org) 
was used to search the subcellular location of each significantly 
changed protein. 

Ingenuity Pathway Analysis 

To identify potential radiation-perturbed molecular pathways in 
the skin secretome, Ingenuity pathway analysis (IPA, www. 
ingenuity.com) was applied to relate the differentially expressed 
proteins to each other based on their interaction and function 
against a high-quality expert-curated knowledge database derived 
from the literature. In brief, all proteins that were statistically 
significandy expressed along with their fold change (log2 ratio) and 
t-test values with respect to 0 cGy group were uploaded into the 
IPA software, which extracted the overlapping networks among 
the candidate proteins from the curated database. Comparison 
analyses were made on both Functions and Canonical Pathways 
levels between different radiation doses, along with a score 
representing the log probability of a particular network being 
found by random chance. The p values were calculated using 
right-tailed Fisher's exact tests. 

Results and Discussion 

Overview of sample preparation, data acquisition and 
data processing 

Identification of secreted proteins from culturing medium 
presents significant challenge due to the composition of condi- 
tioned medium. SDS-PAGE gel analysis showed that significant 
amount of albumin existed in the medium, therefore we used 
immunodepletion to remove the highly abundant medium- 
composing proteins in order to broaden the coverage of the 
proteome secreted by the skin. Depletion of highly abundant 
proteins using avian polyclonal IgY antibody columns has 
demonstrated excellent reproducibility. [26] Therefore, all of the 
medium samples in this study were depleted to enrich the low 
abundant secreted proteins using an IgY- 14 column. 

To overcome the under-sampling issue and increase throughput 
in LC-MS based global proteomics analysis, we applied AMT tag 
strategy for the label free quantitation of each individual 
sample. [17] An on-line 2D LC-MS/MS system was used to 
fractionate and analyze the peptide digests obtained from 4 
samples pooled from each radiation dose group. The peptides 
identified from these LC-MS/MS analyses served as the basis for 
building the skin secretome AMT peptide database. The final 
AMT tag database contained 10,149 peptides from 1547 human 
proteins available for matching to the high resolution LC-MS 
datasets acquired from the individual irradiated media samples. 

After matching the individual sample LC-MS dataset with the 
AMT tag database, 5171 unique peptides were identified and 
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Figure 1. Overview of the A) subcellular location, B) molecular size distribution and C) molecular and cellular functions of the 135 
proteins that showed significant variation between irradiated and sham skin tissues. Others: subcellular location known, but not being 
included in the main categories; Unknown: subcellular location unknown. Protein may have more than one function. 
doi:10.1371/journal.pone.0092332.g001 



quantified from the 48 individual samples. After quality control 
processing to remove outlier LC-MS runs, redundant peptides, 
and peptides with missing data, [23,24] 3386 peptides were 
retained for further downstream data analysis, which include a 
proteomic normalization strategy to select peptide and data scaling 
factor that introduces the least amount of bias into the dataset,[25] 
and comparative statistical analysis was done on the log 
transformed abundance values for each retained peptide to 
calculate the statistical trend, as detailed in the Methods Section. 
Protein roll-up was based on the down-selection of peptides that 
have a common statistical trend followed by a standard average to 
a reference peptide, after filtering the peptides that were 
redundant, had low data content, or were outside the dominant 
significance pattern calculated on peptide level. Trends were based 
on the statistical significance indicator variables devised from 
statistical tests performed on the peptides. The final protein 
quantification was based upon a standard roll-up to a reference 
peptide, where the reference peptide was defined as the peptide 
with the least amount of missing data. In total, 880 human 
proteins were rolled up from 2,229 peptides, and differences in 
protein average abundance between different radiation conditions 
were assessed using a Dunnett adjusted 2-sided t-test to compare 
control group with the radiated groups of 3, 10 and 200 cGy. 
Dunnett p-values less than 0.05 were deemed statistically 
significant. 

Statistical analysis showing dose dependent perturbation 
of secretome by ionizing radiation 

Overall, 135 proteins showed a statistical significant difference 
between the sham (0 cGy) and any of the irradiated groups (3, 10 
or 200 cGy) on the basis of Dunnett adjusted t-test. Thirty three 
proteins with the largest fold change (>1.5 fold for up- and <2.0 
for downregulation) are listed in Table 1 , while the details of all 
proteins are in the Supporting Information (Table SI). To 
facilitate categorization of these 135 proteins, the Uniprot protein 
database was searched to identify their subcellular locations with 
relatively broad definition of subcellular location terms. As shown 
in Figure 1A, majority of these proteins were either originated 
from the cytoplasm (42 %) or direcdy secreted into the extracellular 
matrix (30%), suggesting the quality of our sample preparation 
and analysis process, as well as the potential effects of ionizing 
radiation to the cellular processes. The size distribution of these 



proteins is shown in Figure IB. Functional analysis using IPA 
showed that these proteins were related with protein binding, cell 
signalling, cell movement, epidermis development, endoplasmic 
reticulum stress response, inflammation, free radical scavenging, 
etc. The number of proteins represented by major molecular and 
cellular functions is shown in Figure 1C. 

Among these 135 proteins, 97 proteins showed the trend of 
downregulation with increasing dose of X-rays, irrespective of 
whether they were more or less abundant when compared with the 
sham group (Figure 2A). A few families of proteins dominated this 
category, such as 14-3-3 family, keratin, collagen, alpha-actinin 
family, heat shock protein 70 family, etc. Although for 87 of these 
proteins their changes were only statistically significant (P<0.05) 
for the 200 cGy irradiated group when compared with the sham, 
there were 4 proteins, i.e. putative ubiquitin-conjugating enzyme 
E2 N-like (UE2NL), lactotransferrin (TRFL), secretogranin-3 
(SCG3) and ras-related protein Rab-llA (RB11A) all showed 
significant downregulation for both the 10 and 200 cGy irradiated 
groups (Figure 3A). There were two proteins, endoplasmic 
reticulum resident protein 29 (ERP29) and fibronectin (FINC) 
both upregulated significantly for the 3 and 10 cGy irradiated 
groups (Figure 3A), but the general trend was still going down with 
increasing dose. 

On the other hand, there were 9 proteins that showed the trend 
of upregulation with increasing radiation dose (Figure 2B), 
although for 7 of them the level of upregulation was only 
statistically significant upon exposure to 200 cGy of irradiation. 
Two proteins, carboxypeptidase E (CBPE) and ubiquitin carboxyl- 
terminal hydrolase isozyme Ll(UCHLl) appeared to be sensitive 
to ionizing radiation, but relatively independent to doses received, 
as relatively little change was occurred between different dose 
groups. Conversely, proteasome activator complex subunit 2 
(PSME2) protein was shown to be very sensitive to the dose of 
radiation, as rapid upregulation of this protein was observed when 
radiation doses were increased from 3, to 10 and 200 cGy 
(Figure 3B). 

Unlike most of the proteins that have up or down changes 
associated with increasing radiation dose, there were 21 and 8 
proteins being observed to have the largest change in the 10 cGy 
irradiated group, with the former one having the highest and the 
latter one the lowest level among all three radiated doses, 
respectively (Figure 2C&D). Although irregular, this pattern of 
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Figure 2. Trends associated with protein abundance profiles and doses of ionizing radiation. Protein abundance fold change at each 
radiation dose was normalized to that of sham (0 cGy) using log2 ratio and plotted as the y axis. Each line represents profile of one protein with 
increasing doses of ionizing radiation along the x axis. A) 97 proteins with their abundances decreased; B) 9 proteins with their abundances increased; 
C) 21 proteins with their abundances highest at 1 0 cGy; D) 8 proteins with their abundances lowest at 1 0 cGy. Detailed statistical analysis values, such 
as fold change and p values from Dunnett adjusted 2-sided t-test are listed in Supporting Information Table SI. 
doi:1 0.1 371 /journal.pone.0092332.g002 



change has been observed in the metabolomic profiling of the 
same 3D human skin tissue system exposed to X-ray irradiation as 
well. [7] 

Implications for interpreting the effects of exposure to 
ionizing radiation 

The secretome is defined as the ensemble of proteins secreted 
into the extracellular space by a cell, tissue, organ, or organism at 



any given time and conditions through known and unknown 
secretory mechanisms involving constitutive and regulated secre- 
tory organelles. This term can be broadly defined to include the 
proteins of extracellular matrix (ECM) as well as additional 
proteins shed from the cell surface. Exosomes, due to their small 
size (< 100 nm), are technically challenging to be separated from 
culturing medium. Therefore, the secretome identified in our 
study may contain proteins released from exosomes as well. In 
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Figure 3. Representative proteins that showed A) consistent decrease and B) consistent increase with increasing radiation dose. The 

fold change of protein abundance at each dose was represented by the log2 ratio with respect to the value of sham group. * represents statistically 
significant (P<0.05, Dunnett adjusted 2 sided t-test). Detailed statistical analysis values, such as fold change and p values are listed in Supporting 
Information Table SI. 
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addition to being masked by serum proteins (i.e. fetal bovine 
serum) that normally present in the culture media, for which we 
have applied immunodepletion for their removal, secretome of 
skin tissue could also be contaminated by non-secreted proteins 
released following cell lysis and death. In this respect, the levels of 
two major cytosolic proteins - beta-actin (ACTB_human) and 
beta-tubulin (TBB5_human) have been used by others to monitor 
cell autolysis. [12] Beta-tubulin was not detected in our study. 
Although beta-actin was detected, it didn't appear to be abundant 
on the basis of the abundances of its constitutive peptides; in 
addition, our results did not show any significant changes of this 
protein within the medium samples collected from any of the 
irradiated groups, which demonstrated that the proteomic profile 
changes that we observed were reflective of the true skin 
secretome, not related to ionizing radiation-induced cell death/ 
autolysis. 

Certain forms of soluble keratin were reported to downregulate 
with high doses of X-ray irradiation, [6,2 7] our results concurred 
their findings with the identification of four different forms of 
keratin - K1C10_HUMAN, K2 2 E_HUM AN, K2C 1B_HUMAN, 
K2C1_HUMAN all showed downregulation with respect to the 
sham controls, in particular the downregulation was statistically 
significant at the 200 cGy high dose (see Supporting Information 
Table SI). The 14-3-3 family is comprised of highly conserved 
proteins that are functionally important in the maintenance of 
homeostasis. They directiy interact with various intracellular 
proteins and regulate a large spectrum of both general and 
specialized signalling pathways. [28] Their involvement with the 
cell cycle, association with proto-oncogenes and oncogenes, and 
abnormal expression in various tumors has linked this family of 
proteins to the etiology of human cancer. [29] In this respect, 
tumor suppressor protein 14-3-3-cr was reported to downregulate 
in a radioresistant nasopharyngeal cancer cell line when compared 
with regular nasopharyngeal cancer cell line, [30] and the same 
trend was observed for non-small cell lung cancer cells in 
comparison with normal human bronchial epithelial cells. [31] In 
addition, secreted 14-3-3 zeta has been demonstrated to have 
potential as candidate biomarkers for aggressive thyroid carcino- 
mas. [32] The seven 14-3-3 proteins (see Supporting Information 
Table SI) observed in this work all showed downregulation for all 
of the irradiated samples, in particular for the 200 cGy dose when 
compared with the sham irradiation, which showed that ionizing 
radiation suppress expression of many proteins within the 14-3-3 
family, this may lead to perturbed homeostasis in tumor 
suppression and other important biological functions. 

High doses of ionizing radiation have suppressive effects to 
many biological processes. In this respect, it was reported that 
more downregulated proteins than upregulated proteins were 
being observed after tumor cells exposed to 1 0 Gy of X-rays, as a 
result of cell cycle suspension, suppression of cell proliferation and 
growth, and reduced metabolic processes from the effects of high 
doses of ionizing radiation. [33] Similarly, although measured at 
mRNA level, significant number (244 genes; 33%) of down 
regulated genes was observed when mammary gland tissues were 
exposed to 2 Gy of ionizing radiation; in particular, a number of 
genes with tumor suppressor function persistently downregulated 
in response to radiation exposure. [34] Consistent with these 
reports, majority (102 of 135, Supporting Information Table SI) of 
the significantly changed proteins observed in this study showed 
statistically significant downregulation for the 200 cGy high dose. 
Ingenuity pathway analysis results showed that most of these 



proteins were associated with organization of extracellular matrix, 
cell signalling, cell migration, cell proliferation and apoptosis. 

On the other hand, despite the general downward trend line of 
protein expression with the increase of radiation dose as observed 
in Figure 2A, many of the proteins appeared to be upregulated in 
comparison with sham controls at 3 and 1 0 cGy, although changes 
at these two low doses were largely not statistically significant 
compared to sham controls. This implies that low doses of ionizing 
radiation likely have stimulating effects on the secretome, which is 
in agreement with proteomic studies of endothelial cell exposure to 
low dose gamma radiation (<50 cGy). [35,36] 

Conclusion 

To our knowledge, this work represents the first global 
proteomic profiling and quantitative comparison of the in vitro 
skin tissue secretome after exposure to ionizing radiation. While it 
is consistent with the previous studies that the majority of the 
secreted proteins were downregulated as a result of high dose 
(2 Gy) ionizing radiation-induced stress and damage, due to 
decreased metabolic processes and suppressed cell cycle and 
growth; many of these secreted proteins at low doses were 
upregulated compared with sham controls. Since it is unclear 
whether this response is specific to ionizing radiation or represents 
a general response to injury, the proteins that displayed unique 
patterns of regulation, or being sensitive to either ionizing 
radiation or dose level, warrant further study for their roles in 
adaptive responses of ionizing radiation, and their robustness as 
biodosimetry markers when combined with other agents that may 
induce similar damaging effects. [37,38] 

Data Deposition:. All data related to this work were 
deposited in PeptideAdas, as in http://www.peptideatlas.org/ 
PASS/PASS00404 and http://www.peptideadas.org/PASS/ 
PASS00405. 

Supporting Information 

Table SI Statistically significant proteins secreted into 
the medium by a full thickness skin tissue model at 
different doses (3, 10 and 200 cGy) of X-rays when 
compared with sham controls (0 cGy). Raw data of protein 
abundance were log2 transformed, as represented by the average 
of each irradiated group. Fold changes between treatment and 
sham groups were represented as group differences. Dunnet 
adjusted T-test was used to assess the statistical significance of the 
change. P values <0.05 were deemed statistically significant. 
Dunnet trend: 0, no significant change; 1, upregulation with 
statistical significance; -1, downregulation with statistical signifi- 
cance. Trend: general trend of the average abundance of proteins 
with increasing radiation dose; down, downregulation with 
increasing dose; up, upregulation with increasing dose; A , 10 cGy 
group had the highest abundance; V, 10 cGy group had the lowest 
abundance. Protein name and their subcellular location were 
derived from Uniprot. 
(XLSX) 

Author Contributions 

Conceived and designed the experiments: QZ. Performed the experiments: 
AAS RJM ZH. Analyzed the data: QZ MM BW MEM. Contributed 
reagents/materials/analysis tools: WQRDS WFM. Wrote the paper: QZ. 



PLOS ONE | www.plosone.org 



7 



March 2014 | Volume 9 | Issue 3 | e92332 



Dose Effects of Radiation on Skin Secretome 



References 

1. Royal HD (2008) Effects of low level radiation-what's new? Semin Nucl Med 38: 
392-402. 

2. Hendee WR, O'Connor MK (2012) Radiation risks of medical imaging: 
separating fact from fantasy. Radiology 264: 312-321. 

3. BIER V (2006) Phase 2, Health risks from exposure to low levels of ionizing 
radiation. 

4. Dauer LT, Brooks AL, Hoel DG, Morgan WF, Stram D, et al. (2010) Review 
and evaluation of updated research on the health effects associated with low-dose 
ionising radiation. Radiat Prot Dosimetry 140: 103—136. 

5. Brenner DJ, Doll R, Goodhead DT, Hall EJ, Land GE, et al. (2003) Cancer risks 
attributable to low doses of ionizing radiation: assessing what we really know. 
Proc Natl Acad Sci U S A 100: 13761-13766. 

6. Yang F, Waters KM, Webb- Robertson BJ, Sowa MB, von Neubeck C, ct al. 
(2012) Quantitative phosphoproteomies identifies filaggrin and other targets of 
ionizing radiation in a human skin model. Exp Dermatol 21: 352—357. 

7. Hu ZP, Kim YM, Sowa MB, Robinson RJ, Gao X, et al. (2012) Metabolomic 
response of human skin tissue to low dose ionizing radiation. Mol Biosyst 8: 
1979-1986. 

8. Snijders AM, Marchetti F, Bhatnagar S, Duru N, Han J, et al. (2012) Genetic 
differences in transcript responses to low-dose ionizing radiation identify tissue 
functions associated with breast cancer susceptibility. PLoS ONE 7: e45394. 

9. Knops K, Boldt S, Wolkenhaucr O, Kriehuber R (2012) Gene expression in low- 
and high-dose -irradiated human peripheral blood lymphocytes: possible 
applications for biodosimetry. Radiat Res 178: 304-312. 

10. Chaudhry MA, Omaruddin RA, Krcger B, de Toledo SM, Azzam EI (2012) 
Micro RNA responses to chronic or acute exposures to low dose ionizing 
radiation. Mol Biol Rep 39: 7549-7558. 

11. Albrccht H, Durbin-Johnson B, Yunis R, Kalanetra KM, Wu S, ct al. (2012) 
Transcriptional response of ex vivo human skin to ionizing radiation: 
comparison between low- and high-dose effects. Radiat Res 177: 69-83. 

12- Makridakis M, Vlahou A (2010) Secretome proteomics for discovery of cancer 
biomarkers. J Proteomics 73: 2291-2305. 

13. Morgan WF, Sowa MB (2007) Non-targeted bystander effects induced by 
ionizing radiation. Mutat Res 616: 159-164. 

14. Varnum SM, Springer DL, Chaffee ME, Lien KA, Webb-Robertson BJ, et al. 
(2012) The effects of low-dose irradiation on inflammatory response proteins in a 
3D reconstituted human skin tissue model. Radiat Res 178: 591-599. 

15. Mothcrsill C, Antonelli F, DahleJ, Dini V, Hegyesi H, ct al. (2012) A laboratory 
inter-comparison of the importance of serum serotonin levels in the 
measurement of a range of radiation- induced bystander effects: overview of 
study and results presentation. Int J Radiat Biol 88: 763-769. 

16. Pichowski PD, Pctyuk VA, Orton DJ, Xie F, Moore RJ, ct al. (2013) Sources of 
technical variability in quantitative LC-MS proteomics: human brain tissue 
sample analysis. J Proteomc Res 12: 2128—2137. 

17. Zimmcr JS, Monroe ME, Qian WJ, Smith RD (2006) Advances in proteomics 
data analysis and display using an accurate mass and time tag approach. Mass 
Spectrom Rev 25: 450-482. 

18. Mayampurath AMJaitly N, Purvine SO, Monroe ME, Auberry KJ, et al. (2008) 
DeconMSn: a software tool for accurate parent ion monoisotopic mass 
determination for tandem mass spectra. Bioinformatics 24: 1021-1023. 

19. Kim S, Gupta N, Pcvzner PA (2008) Spectral probabilities and generating 
functions of tandem mass spectra: a strike against decoy databases. J Proteomc 
Res 7: 3354-3363. 

20. Petritis K, Kangas LJ, Yan B, Monroe ME, Strittmatter EF, et al. (2006) 
Improved peptide elution time prediction for rcverscd-phase liquid chromatog- 
raphy-MS by incorporating peptide sequence information. Anal Chem 78: 
5026-5039. 

21. Jaitly N, Mayampurath A, Littlefield K, Adkins JN, Anderson GA, et al. (2009) 
Decon2LS: An open-sourec software package for automated processing and 



visualization of high resolution mass spectrometry data. BMC Bioinformatics 10: 

87. 

22. Monroe ME, Tolic N, Jaitly N, ShawJL, Adkins JN, et al. (2007) VIPER: an 
advanced software package to support high-throughput LC-MS peptide 
identification. Bioinformatics 23: 2021-2023. 

23. Wcbb-Robcrtson BJ, McCue LA, Waters KM, Matzke MM, Jacobs JM, et al. 
(20 1 0) Combined statistical analyses of peptide intensities and peptide 
occurrences improves identification of significant peptides from MS-based 
proteomics data. J Proteomc Res 9: 5748-5756. 

24. Matzke MM, Waters KM, Metz TO, Jacobs JM, Sims AG, et al. (2011) 
Improved quality control processing of pcptidc-centric LC-MS proteomics data. 
Bioinformatics. 

25. Wcbb-Robcrtson BJ, Matzke MM, Jacobs JM, Pounds JG, Waters KM (201 1) A 
statistical selection strategy for normalization procedures in LC-MS proteomics 
experiments through datasct-dependent ranking of normalization sealing factors. 
Proteomics 11: 4736-1741. 

26. Shi T, ZhouJY, Gritsenko MA, Hossain M, Camp DG, 2nd, et al. (2012) IgY14 
and SupcrMix immunoaffinity separations coupled with liquid chromatography- 
mass spectrometry for human plasma proteomics biomarkcr discovery. Methods 
56: 246-253. 

27. Nanashima N, Ito K, Ishikawa T, Nakano M, Nakamura T (2012) Damage of 
hair follicle stem cells and alteration of keratin expression in external radiation- 
induced acute alopecia. Int J Mol Med 30: 579-584. 

28. Mhawech P (2005) 14-3-3 proteins-an update. Cell Res 15: 228-236. 

29. Wang Z, Ncsland JM, Suo Z, Trope CG, Holm R (201 1) The prognostic value 
of 14-3-3 isoforms in vulvar squamous cell carcinoma cases: 14-3-3beta and 
epsilon arc independent prognostic factors for these tumors. PLoS ONE 6: 
c24843. 

30. Feng XP, Yi H, Li MY, Li XH, Yi B, ct al. (2010) Identification of biomarkers 
for predicting nasopharyngeal carcinoma response to radiotherapy by 
proteomics. Cancer Res 70: 3450-3462. 

31. Liu Y, Chen ZhangJT (2004) Tumor suppressor gene 14-3-3sigma is down- 
regulated whereas the proto-oncogene translation elongation factor ldelta is up- 
rcgulated in non-small cell lung cancers as identified by proteomic profiling. 
J Proteomc Res 3: 728-735. 

32. Chakcr S, Kashat L, Voisin S, KaurJ, Kak I, ct al. (2013) Secretome proteins as 
candidate biomarkers for aggressive thyroid carcinomas. Proteomics 13: 771- 
787. 

33. Wang F, Bing Z, Zhang Y, Ao B, Zhang S, ct al. (2013) Quantitative proteomic 
analysis for radiation-induced cell cycle suspension in 92-1 melanoma cell line. 
J Radiat Res 54: 649-662. 

34. Datta K, Hydukc DR, Suman S, Moon BH, Johnson MD, ct al. (2012) Exposure 
to ionizing radiation induced persistent gene expression changes in mouse 
mammary gland. Radiat Oncol 7: 205. 

35. Barjaktarovic Z, Anastasov N, Azimzadeh O, Sriharshan A, Sarioglu H, et al. 
(2013) Integrative proteomic and micro RNA analysis of primary human 
coronary artery endothelial cells exposed to low-dose gamma radiation. Radiat 
Environ Biophys 52: 87-98. 

36. Pludcr F, Barjaktarovic Z, Azimzadeh O, Mortl S, Kramer A, et al. (201 1) Low- 
dose irradiation causes rapid alterations to the proteome of the human 
endothelial cell line EA.hy926. Radiat Environ Biophys 50: 155-166. 

37. Amundson SA, Bittner M, Meltzer P, Trent J, Fornacc AJ, Jr. (2001) Biological 
indicators for the identification of ionizing radiation exposure in humans. Expert 
Rev Mol Diagn 1: 211-219. 

38. Laiakis EC, Hyduke DR, Fornacc AJ (2012) Comparison of mouse urinary 
metabolic profiles after exposure to the inflammatory stressors gamma radiation 
and lipopolysaccharide. Radiat Res 177: 187-199. 



PLOS ONE | www.plosone.org 



8 



March 2014 | Volume 9 | Issue 3 | e92332 



